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Abstract 

The electronic and magnetic structures of ScFe2 and of its dihydride ScFe2H2 are self-consistently 
calculated within the density functional theory (DFT) using the all electron augmented spherical 
wave (ASW) method with the local spin density approximation (LSDA) for treating effects of 
exchange and correlation. The results of the enhancement of the magnetization upon hydrogen 
^ I insertion are assessed within an analysis of the chemical bonding properties from which we suggest 

'^ ■ that both hydrogen bond with iron and cell expansion effects play a role in the change of the 

O . magnitude of magnetization. In agreement with average experimental findings for both the inter- 

metallic system and its dihydride, the calculated Fermi contact terms Hpc of the ^"^Fe Mossbauer 
^ ■ spectroscopy for hyperfine field, at the two iron sites, exhibit an original inversion for the order of 

Cn . magnitudes upon hydriding. 
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I. INTRODUCTION 

Binary alloys belonging to the Laves family AB2 (A= rare earth or actinide, B = transition 
metal) crystallize either in a face centered cubic fee (C15) lattice or in a hexagonal lattice 
(C14); further, dihexagonal (C36) minority structure exists.-*^ Within this family, ScFe2 
exhibits polymorphism and can exist in the three above crystalline states.- One of the basic 
aspects of the electronic properties of the pure AB2 alloy systems is the identification of the 
origin of the magnetism which can be either due to the transition metal or induced by the 
A metal, depending on the chemical nature of the involved species.-i^ Among others,-i^ Smit 
and Buschow have studied the synthesis of the ScFe2 compound^, reporting experimental 
measurements for the average magnetic moment of iron and the effective hyperfine field 
Heff. On the other hand, the interaction of these intermetallic phases with hydrogen was 
investigated in a number of works. -i^i^^iii Besides the large potential applications of Laves 
phases hydrides in the field of solid state storage of hydrogen for energetics,— ^^^ there is 
a basic interest in studying the magnetic structure and the electronic properties due to 
H insertion. At this level of investigation concerning ScFe2, ^'^Fe Mossbauer spectroscopic 
works show an increase for the average magnitudes of the magnetization and the hyperfine 
field for iron upon hydriding without specifically assigning a role for each one of the two 
iron sites. ^ This leads to suggest an interplay between magnetovolume and chemical effects 
brought by the cell expansion when hydrogen is inserted. This original feature is addressed 
in this work. Further, a detailed atom-resolved study of the magnetism is provided and the 
nature of the non-rigid-band behavior within ScFe2 and its dihydride are assessed. 

II. CRYSTAL STRUCTURES 

The ScFe2 intermetallic system is experimentally stable in the hexagonal C14-type, 
PG^/mmc space group. Laves structure.- For this alloy, Sc atoms are located in 4f sites 
at (1/3, 2/3, 0.0661). As for Fe atoms, there are two crystallographic nonequivalent sites: 
Fel(2a) at (0, 0, 0) and Fe2(6h) at (0.8357, 1.6714, 1/4). Fe2, Sc and Fel atoms have an 
occupancy ratio of 3 : 2 : 1. All sites mentioned in this report are in Wyckoff notation; 
numerical values are given in Refs.— li^ii^. Smit and Buschow^ charged an alloy sample with 
hydrogen gas, the composition of the formed hydride was found to correspond to ~ 1.92 



atoms of H per formula unit (fu). Furtheron, this hydride will be referred to as ScFe2H2. A 
sketch of the structure containing hydrogen is presented in Fig. [H 




FIG. 1: The hexagonal crystal structure of ScFe2H2 ('C14' Laves phase, space group PQ^/mrac). 
Sc (light grey), Fel (grey), Fe2 (black) and H (small light grey) are drawn in decreasing sphere 



sizes. 



The choices for hydrogen insertion sites were done based on the neutron diffraction studies 
of Didisheim on the deuterated C14 compound ZrMn2.J^ Those studies show that hydrogen 
atoms in ternary hydrides, based on the C14-type, occupy interstitial tetrahedral-like A2B2 
holes in positions such as: 241, 12k, 6hl and 6h2. Smit, Donkersloot and Buschow^^, have 
hydrided ScFe2 into ScFe2H2 based on that same study, suggesting that each type of these 
holes is partially occupied by H atoms. In the present work, calculations were done with the 
assumption that hydrogens are located in sites such as: Hl(6hl) at (0.463, 0.926, 1/4) and 
H2(6h2) at (0.202, 0.404, 1/4). There are 12 of these tetrahedral holes per fu. Nevertheless, 



H atoms will occupy only 8 of these interstices with respect to the stoichiometry of Sc4Fe8H8. 
The occupancy ratios of HI and H2 atoms are the same as those of Fe2 and Fel in ScFe2 
respectively. Moreover, 12k and 241 interstices were tested. The performed calculations 
showed that they cannot reproduce magnitudes for the magnetic moments that are coherent 
with the expreriment. Values of the average magnetic moment for iron such as 1.865 and 
2.065 fiB, for 12k and 241 respectively, fail to compare with the experimental value of 2.23 
Hb-— These two A2B2 holes were discarded in favor of 6hl and 6h2 intersitices which give 
more suitable values for the magnetic moments that agrees with experiment; these will be 
given later within the text. One can notice upon examining Fig. [T]that Fel atoms are little 
bonded with H atoms. As a matter of fact, A2B2 holes are formed by two 3d transition 
metal atoms (Fel and Fe2) and two atoms of the strongly hydrogen attracting component 
(Sc), which explains H atoms bonds with Sc. As for H-Fe2 bonds, table [TTl reports shorter 
Fe2-H distances compared to Sc-H one, which favors the H intake into A2B2 sites formed 
by Fe2 with respect to those formed by Fel. Keeping in mind the relatively large difference 
in atomic volume between Sc and Fe, the choice of the H insertion sites is in agreement with 
the Westlake criterion (— and therein cited Refs.) that imposes a minimum interstitial hole 
size of 0.40 A. Too small H-H distances {(Ih^h > 2.lA) were avoided in order to respect 
the Switendick criterion.-- Table [T] provides the crystal data of the two systems discussed in 
Ref.-. Note that there is a mismatch of the c/a ratio magnitudes for both experimentally 
prepared ScFe2 and ScFe2H2 as stated in Ref.-, i.e., c/a = 1.636 for ScFe2 and c/a = 1.611 
for ScFe2H2. Thus an anisotropic evolution accompanies the formation of the dihydride. 

In this work hydrogen insertion effects within ScFe2 are examined within three com- 
plementary approaches relevant (i) to the influence of volume expansion on the magnetic 
properties within the structure, (ii) to the chemical role of hydrogen and its influence on 
the changes of the magnetization for the dihydride lattice and (iii) to the crystal anisotropy 
occurring upon hydriding the alloy system. 





ScFe2 


Expanded 


ScFe2 


ScFeaHs 




c/o=1.636 


c/a=1.636 


c/a=1.611 


c/a=1.611 


a (A) 


4.963 


5.250 


5.277 


5.277 


c(A) 


8.122 


8.592 


8.504 


8.504 


Volume 


43.321 


51.276 


51.276 


51.276 


dsc-Fel 


2.915 


3.084 
3.110 


3.100 


3.100 


dSc-Fe2 


2.899 


3.052 


3.068 


3.068 




2.952 


3.084 


3.100 


3.100 


dFel-Fe2 


2.470 


2.592 
2.698 


2.603 


2.603 


dm-Sc 








1.963 


dHl-Fe2 








1.724 


dH2-Sc 








1.973 


dH2-Fe2 








1.698 



TABLE I: ScFe2 model systems: Crystal data for ScFe2 and ScFe2H2^ are presented in columns 1 
and 4 respectively. Volume is given in {A? / fu ). For more details see text. 

III. COMPUTATIONAL FRAMEWORK 

A. Electronic and magnetic properties 

Among the self-consistent methods built within the density functional theory (DFT) 
(see for instance^'^'^) we use the augmented spherical wave (ASW) method^"^ which is 
an "all-electrons" method (non frozen core). The ASW method has proven its efficiency in 
similar topics of intermetallic systems and their hydrides.— i^^ The analysis of the calculation 
results allows assigning a role to each atomic constituent in the magnetism and in the 
chemical bonding. The effects of exchange and correlation were treated within a local spin 
density approximation (LSDA) scheme.— All valence electrons were treated as band states. 
In the minimal ASW basis set,— we chose the outermost shells to represent the valence 
states and the matrix elements were constructed using partial waves up to Imax. + 1 = 3 



for Sc and Fe and Imax. + 1 = 2 for H. The completeness of the valence basis set was 
checked for charge convergence meaning that charge residues are < 0.1 for Imax. + 1- The 
self-consistent field calculations were run to a convergence of AQ = 10^^ for the charge 
density^^ and the accuracy of the method is in the range of about 10"'^ eV regarding energy 
differences. Besides its construction within the DFT, the ASW method is based on the 
atomic sphere approximation (ASA) which assumes overlapping spheres centered on the 
atomic sites within which the potential has a spherical symmetry -central potential-. The 
volume of the spheres has to be equal to the cell volume because the wave equation is solved 
only in the spheres. This is unproblematic for closely packed structures like metals and 
intermetallics such as ScFe2 itself. But for less compacked structures such as that of ScFe2H2 
studied here, additional augmentation spheres, called empty spheres (ES) are introduced to 
represent the interstitial space without loss of crystal symmetry and to avoid an otherwise 
too large overlap between the actual atomic spheres. ES are "pseudo atoms" with zero 
atomic number. They receive charges from the neighboring atomic species and allow for 
possible covalency effects within the lattice. Within the ASW method, the sphere geometry 
optimization (SGO)^ algorithm is used to generate ES without symmetry breaking. For 
ScFe2H2 one type of ES has been added. In total, 12 ES have been inserted into the Sc4Fe8H8 
structure. Besides obtaining the electronic band structure with site projected density of 
states (PDOS) for total spins (non spin polarized NSP configuration), as well as the spin- 
resolved PDOS, the calculations allow discussing quantities such as the magnetic moments 
and their sign, the magnitude of the exchange splitting, as well as the spin densities at the 
core due to the ns polarization thanks to the d magnetic moment. These lead to the Fermi 
contact term of the hyperfine field Hpc which constitutes the major part of the hyperfine 
field obtained by ^^Fe Mossbauer spectroscopy. 

The calculations are started by assuming a non-magnetic configuration meaning a spin 
degeneracy for all valence states and equal spin occupations. Such a configuration should not 
be confused with a paramagnet, which could be simulated either by a supercell calculation 
with random spin orientations or by calling for disordered local moment approaches based 
on the coherent potential CPA approximation^^ or the LDA+DMFT scheme.— Subsequent 
spin-polarized calculations with different initial spin populations can lead at self-consistency 
either to finite or zero local moments within an implicit long-range ferromagnetic order. 



B. Chemical bonding properties 

The interactions within the alloy lattice with inserted hydrogen can be described in the 
framework of chemical bonding. Some elaborate tools exist allowing to obtain information 
about the nature of such interactions between atomic constituents as well as the respective 
quantum states involved. This can be provided using overlap population Sij (OP) leading 
to the so-called crystal orbital overlap population (COOP)^^ or alternatively introducing 
the Hamiltonian based population {Hij) with the crystal orbital Hamiltonian population 
(COHP).— Both approaches provide a qualitative description of the chemical interactions 
between two atomic species by assigning a bonding, non-bonding or anti-bonding character. 
A slight refinement of the COHP was recently proposed in form of the "covalent bond en- 
ergy" ECOV which combines both COHP and COOP so as to make the resulting quantity 
independent of the choice of the zero of potential.— Our experience with both COOP— and 
ECOV^'^ shows that they give similar general trends although COOP exaggerate the mag- 
nitude of anti-bonding states. The ECOV criterion implemented within the ASW method 
is used here for the description of the chemical bond. 

IV. ScFea VERSUS VOLUME EFFECTS 

Calculations were performed at the experimental volume of ScFe2.- In order to evaluate 
the magnetovolume effects, the computed results for the expanded hydrogen free ScFe2 at 
the same lattice constants of the dihydride were addressed. As a matter of fact, such effects 
can be important in these intermetallic systems in as far as the onset of the magnetic moment 
is due to interband spin polarization, i.e., it is mediated by the electron gas in a collective 
electrons approach. This is opposite to other systems, such as insulating oxydes where the 
magnetization is of intraband character, and hence, less affected by volume changes such as 
those induced by pressure (negative or positive).— 

A. NSP calculations 

At self consistent convergence a progressive increase of BZ integration up to the value of 
576 k-points, i.e., 16 k-points in each direction of the irreducible wedge of the hexagonal BZ 
was used. A slight charge transfer of ~ 0.104 electron is seen from Fe2 towards Sc and Pel. 



However its amount is not significant of an ionic beliavior -rarely observed in tlie framework 
of ab initio calculations for such systems—. Therefore it can be argued that the bonding 
is not mainly due to charge transfer but rather imposed by the hybridization between the 
different valence states. It is also important to mention that for all calculations (NSP as 
well as SP) the best evaluation of the radii for the different atomic species was assumed 
resulting in a better ASA overlap. 
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FIG. 2: Non magnetic site projected DOS of ScFe2 



1. Projected density of states PDOS 

The PDOS for the ScFe2 is given in Fig. [2] with respect to occupancy ratios given in 
section II (twice more Fe than Sc). This is applied for all the other PDOS pannels in this 
work. The origin of energies along the x axis is taken with respect to the Fermi energy 
{Ep)] this is equally followed in all other plots. Looking firstly at the general shape of the 
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PDOS one can observe that the Ep level is situated at the peaks of both Fel and Fe2 with 
a predominance in terms of intensity for Fe2 states, i.e., with respect to very low intensity 
scandium states. The similar skylines between the partial PDOS pointing to the mixing 
between Fe2, Fel and Sc states can be seen at the lower part of the valence band (VB), 
with mainly sp like states between —6 and —2.5 eV, as well as towards the top of VB {d 
states). Such mixing will be analyzed later regarding the chemical bonding. Lastly, within 
the conduction band (CB), 3(i(Sc) states are found dominant. This is expected as scandium 
is located at the very beginning of the 3d period, with mainly empty d states. 

2. Analysis of the NSP-PDOS within the Stoner theory 

In as far as Fel, Fe2 and Sc 3d states were treated as band states by our calculations, 
the Stoner theory of band ferromagnetisn>2i can be applied to address the spin polarization. 
The total energy of the spin system results from the exchange and kinetic energies counted 
from a non-magnetic state. Formulating the problem at zero temperature, one can express 
the total energy as E = ^[ J^ J [l — ln{Ep)]. Here I is the Stoner exchange-correlation 
integral which is an atomic quantity that can be derived from spin polarized calculations.— 



n{Ep) is the PDOS value for a given species at the Fermi Level in the non-magnetic state. 
The product \n{EF) from the expression above provides a criterion for the stability of the 
spin system. The change from a non-magnetic configuration towards spin polarization is 
favorable when \n{Ep) > 1. The system then stabilizes through a gain of energy due to 
exchange. From Ref.— , I(Fe) = 0.4624 eV and the computed n{Ep) values for Fel and Fe2 
are ~ 3.047 and 3.286 eV~^ respectively. The Stoner products for Fel and Fe2 are then 
~ 1.409 and ~ 1.519 respectively. This means that the Stoner criterion is satisfied for the 
two iron sites, within ScFe2. 

3. Covalent bond energy ECOV 

The analysis of the chemical bonding is done using the ECOV approach presented in 
section III.B . The corresponding plots are shown in Fig. [SJ 

Along the coordinate axis negative, positive and zero ECOV magnitudes (unitless) point 
to bonding, anti-bonding and non-bonding interactions respectively. Partial ECOV are 
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FIG. 3: Chemical bonding: non magnetic ECOV for ScFe2. 

given for the atomic pair interactions of Sc-Fel, Sc-Fe2 and Fel-Fe2 bonds. The other 
ECOV plots presented in this work equaly describe atomic-pair interactions. The change 
in bonding strength is proportional to the distance magnitudes given in table [H i.e., the 
shortest interatomic distances characterize the strongest interactions. The anti-bonding 
character of the strongest Fel-Fe2 interaction at the top of the VB and at Ep points to 
the instability of the system in the non-magnetic NSP configuration. On the contrary, the 
bonding character of the Sc-Fel interaction up to Ep contributes to the stability of the 
system. The d band electrons crossed by the Fermi level are not all anti-bonding. A part of 
those electrons becomes non-bonding in the neighborhood of Ep, thus participating to the 
onset of the magnetic moment. 
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FIG. 4: Site and spin projected DOS of ScFf^ (a) and the expanded hydrogen free ScFe2 (b) 
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FIG. 5: Spin resolved chemical bonding: | and J, ECOV-spin for Sc-Fel atomic pair interaction 
within ScFe2. 

B. Spin polarized calculations 

As it can be expected for the magnetic configuration of ScFe2, there is an energy stabi- 
hzation of AE = 0.218 eV per fu with respect to the NSP calculations. This agrees with 
the experimental ferromagnetic ground state whereby ScFe2 is considered as the only ferro- 
magnet among C14 stoichiometric transition metals compounds as described in Refs.-'^'^. 

1. Magnetic moments 

Magnetic moments are obtained from the charge difference between | —spin and J, —spin 
of all valence states; their calculated values are listed in table [TTl The computed val- 
ues for both the average magnetic moment of iron and the magnetization per fu are 
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ScFes 


Expanded 


ScFe2 


ScFeaHs 




c/a=1.636 


c/a=1.636 


c/a=1.611 


c/a=1.611 


msc 


-0.482 


-0.676 


-0.682 


-0.347 


mpel 


1.468 


2.229 


2.279 


2.486 


mFe2 


1.560 


2.101 


2.118 


1.967 


< mpe > 


1.514 


2.165 


2.198 


2.226 


M 


2.600 


3.591 


3.635 


3.791 


H'^'c^\Fel) 


-157 


-240 


-246 


-203 


H*°if\Fe2) 


-162 


-229 


-226 


-232 


H^^^%Fel) 


-159 


-252 


-255 


-272 


Hf^%Fe2) 


-164 


-232 


-231 


-214 


Erel 


24.733 


26.661 


26.671 


0.000 



TABLE II: Magnetic results for ScFe2 calculated in this work. Magnetic moments and total values 
are given in ^b- The core part H'j^^^ of the fermi contact term of Heff a-s well as H^p^^ which 
is the sum of both core and valence contributions have their calculated value listed in kGauss. 
Ej.(,i represents the relative SP-energy per fu towards the most stable value of —89884.16355 eV 
corresponding to ScFe2H2. For more details see text. 

< TTiFe >= 1.514 fiB and M = 2.600 fiB- Different experimental values were obtained by 
magnetic measurements; M = 2.9 yU^^^ and < rripe >= 1.2, 1.37, 1.45 fiB given in Refs.— i^Ii^ 
respectively. Smit and Buschow also reported in Ref.— that < rupe > is equal to 1.14 fiB 
for sub-stoicliiometric ScFei.gg and to 1.34 fiB for over-stoichiometric ScFe2.o5- It can be 
then suggested that the values in this work are within the range of the experimental data. 
Also the scandium carries a negatif magnetic moment of —0.482 fiB- From this ScFe2 is a 
ferrimagnet in its ground state, contrary to some experimental results^*^ which announce it 
as a ferromagnet. The same ordering was observed for yttrium, within YFe2, which carries 
a magnetic moment of —0.50 fiB-— Moreover, the 3(i(Sc) orbital holds a calculated value 
of —0.347 fiB which stands out as the largest contribution within the magnetic moment. 
The second largest contribution is that of the 4p(Sc) orbital with a value of —0.103 fiB- 
On the other hand, one can establish on analyzing the magnetic results for the expanded 
hydrogen free ScFe2 given in table [TTl that ScFe2 gives rise to an increase in the magnetic mo- 
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ments {Am/m) for all the atomic species upon volume expansion. This increase is around 
64% for the average magnetic moment of iron. Experimental values^ are found such as 
Am/m = 0.96 and 0.37 for ScFei.gg and ScFe2.o5 respectively. The fact that the calculated 
value ranges between these two experimental values is related to the stoichiometric ordered 
lattice assumed by the calculations. Furthermore, the values Am,/m = 0.46,0.75 and 0.87 
for Fe2, Sc and Fel respectively show that this increase is not at the same rate for the differ- 
ent species. This can be explained in terms of a reduction in contact between the scandium 
and iron atoms.— The interatomic distances values in tabled confirm this explanation where 
it is found that the distance between the atomic pair Sc-Fe2 is smaller than the Sc-Fel dis- 
tance. This is opposite to the ScFe2 alloy system where d^c-Fei < '^Sc-Fe2 and m,pei < fnFe2- 
Lastly, table [III reports a higher mpe2 magnitude with respect to m^p^i for the intermetallic 
ScFe2. This order of magnitudes is not respected for the expanded hydrogen free ScFe2, 
where m,Fe2 is smaller than m^p^i. The calculations show an increase in the difference be- 
tween I and I —spin for the d states upon volume expansion for both Fel and Fe2. This 
increase is around 85 and 46% for Fel and Fe2 respectively. This major difference explains 
the inversion for the order of magnitudes brought by volume expansion. 

2. Hyperfine field Hpc 

Another significant result extracted from these calculations is the Fermi contact term 
of the hyperfine field (Hpc) (see table [TTl)- The effective magnetic field H^jf acting on a 
nucleus is considered as the sum of four contributions; (i) Hi, the internal field which is the 
magnetic field at the nucleus generated from an externally applied field, (ii) Hpc, the Fermi 
contact term, based on the spin density at the nucleus for the ns quantum states caused by 
the polarization of the s electrons by the d moments, (iii) Horb, which is the field arising 
from the orbital magnetic moment and (iv) Hdip, representing the dipole interaction with 
the surrounding atoms. In a non-relativistic description, Hpc is expressed by the formula: 
Hpc = — ^7Af{($|(0))^ — ($|(0))^}. Where 7Ar is the nuclear gyromagnetic ratio and the 
quantities between brackets are the densities of s electrons at the nucleus (r = 0) for | 
and I —spin respectively. The calculated Hpc values are such as: Hp^\Fel) = —157 
kGauss and H^p^\Fe2) = —162 kGauss. Based on ^^Fe Mossbauer spectroscopy, reported 
experimental values of H^f^ are —167 and —174 kGauss for Fel and Fe2 respectively. 
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Other experimental values of the average hyperfine field were observed to be —170 kGauss 
for Fel and —176 kGauss for Fe2 for ScFei.ge.— Magnitudes of —220 and —203 kGauss, for 
the two crystallographic iron sites 2a and 6h respectively, are found for ScFe2.o5 in Ref.— . 
This is an experimental evidence of the inversion for the order of magnitudes of H^ff in the 
intermetallic ScFe2 due to stoichiometric changes. Such small departures from stoichiometry 
cannot be studied in the scheme of the calculations performed in this work; they require 
other schematic representations such as with the CPA.— The difference between calculated 
and experimental values can be related to different origins relevant to (i) the fact that 
the local spin density approximation cannot treat with sufficient accuracy the polarization 
of core wave functions,^^ (ii) the non-stoichiometry of the experimetally prepared alloys 
and the subsequent disorder within the solid solutions. An explanation for this peculiar 
behavior, i.e. the difference between calculated and measured values, can be found by 
decomposing Hpc into its major contributions namely the one from the core Is, 2s and 
3s electrons H'^q'^ and the one from the valence 4s electrons Hp"/^.— While H'^q'^ usually is 
strictly proportional to the magnetic moment, Hp'q contains large contributions from the 
neighboring atoms. The calculated H^q'' values presented whitin this work are the sum of 
these two parts. Considering only the H'^^'^ contribution, values such as —159 and —164 
kGauss were obtained for Fel and Fe2 respectively. This is a slight enhancement of the 
calculations with respect to experimental measurements. 

3. Projected density of states PDOS 

The PDOS curves for the spin polarized SP configuration of ScFe2 are shown in Fig. 111(a). 
Within the VB two energy regions can be identified, from —7.5 to —5 eV, low intensity 
itinerant s,p states of all constituents are found; this is followed by larger intensity peaks 
mainly due to 3d{Fe) up to and above Ep. Exchange splitting can be seen to mainly affect 
the latter as it is expected from the above analysis of the magnetizations. | —spin states, 
for both Fel and Fe2 at Ep, are concentrated in sharp and narrow PDOS peaks, contrary 
to I —spin states that are found in PDOS minima. Fel and Fe2 peaks within the energy 
range [—1.3, —0.3] eV for f —spin states are similar to those corresponding to [0, 1] eV for 
I —spin states. This shift in spectral weight for | —spin states below the Fermi level and 
for J, —spin states above Ep corresponds to the onset of magnetic moments carried by Fel 
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and Fe2. One can attribute this to a Stoner rigid-band magnetism at first sight. But PDOS 
weights at | and J, —spin populations are not the same. This mismatch between both spin 
populations is mainly due to the 3d{Fe2) states peaks at ~ 1 eV for the | —spin states. 
This implies a non rigid-band shift which rules the magnetism of this system unlike aFe.— 
Magnetism arising in this way is called "covalent magnetism" M Such behavior was formerly 
shown for ZrFe2.— One also notices that the peaks at ~ —0.7 eV for | —spin 3d{Sc) states 
are more intense than those for | —spin states, a magnetic moment is carried by Sc with 
smaller magnitude and opposite direction to those carried by Pel and Fe2. The moment of 
scandium is provided by the covalent Sc-Fel bond, rather than by a rigid energy shift of 
non-magnetic PDOS, whence its negative sign -notice the Sc-Fel overlap around —0.7 eV 
for I —spin PDOS -. On the other hand, the SP-PDOS for the expanded hydrogen free 
ScFe2, given in Fig. HI show through the intense peaks at ~ —0.5 eV, a larger | —spin 
states occupation near Ep in respect with the j —spin states. This feature, brought by 
volume expansion, is responsible of the inversion for the order of magnitudes of magnetic 
moments (see section IV.B.l). Another important feature is the Sc intense peak for J, —spin 
at ~ —0.5 eV which is shifted closer to Ep and is found more intense with respect to the 
scandium peak at ~ —0.7 eV in Fig. Hl^b). This peak is responsible for the increasing of the 
magnetic moment for Sc (see section IV.B.l) which is related to volume expansion. 

4- Spin resolved chemical bonding 

The Sc-Fel interactions for J, and t —spin states curves are given in Fig. [51 The plots 
show that the | —spin bonding is stronger with respect to the | —spin one. This points to 
spin resolved chemical bonding. A peak at ~ —0.7 eV is pointing for the Sc-Fel | —spin 
interactions confirming this bond contribution to the magnetic moment carried by scandium. 

V. ScFea AND HYDROGEN EFFECTS 

For ScFe2H2, computed at the experimental data given in Ref.-, the BZ integration within 
the self-consistent calculations was carried out up to 1024 k -points. Charge transfer within 
ScFe2H2 shows a departure of ~ 0.55 electron from Sc spheres to Pel, Fe2 and empty spheres. 
This slight transfer, not significant of ionic effects, signals a redistribution of the two s 
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electrons of Sc over its three valence basis sets thus providing it with a larger d character 
arising from its mixing with 3(i(Fel) and 3(i(Fe2). The covalent bond between the metal 
species and H can be rather dicussed using the PDOS and ECOV criteria. The magnetic 
configuration of the dihydride is, as expected, energetically stabilized by ^E = 0.436 eV 
and 15.451 eV per fu compared to the NSP calculations performed for ScFe2H2 and the SP 
calculations for ScFe2 respectively. The relative energies per fu (Erei) for the magnetic SP 
configuration, reported in table [TTl describe the difference between the most stable ScFe2H2 
energy and the energies of the other computed systems. The values for Erei show that while 
the volume expansion destabilizes the system, the hydriding restores it to an even more stable 
state. Using an equation giving the binding energy of the hydride system follwoingi^ for 2H 
per fu: Eb = \ [-E'(ScFe2H2) — -E'(ScFe2) — -E'(H2)], a calculation for the binding energy of 
the system due to hydrogen is performed and the resulting value of —0.92 eV.mol"^ is found 
in agreement with values for other systems.— This shows the important role of hydrogen in 
the stabilization of the intermetallic system. 

A. Magnetic moments 

Magnetic moments values derived from spin polarized calculations are given in table [TTl 
The average magnetic moment for iron is of 2.226 /is which agrees with the experimental 
result of 2.23 ^b given in Ref.^. This computed value corresponds to the increase I\mjm = 
0.67, which is an enhancement of 1.3% over the increase brought by the volume expansion. 
Moreover, the respective Am/m values for metal species (with respect to ScFe2) show a 
decrease of 11% for Sc as well as an increase of 35 and 104% for Fe2 and Fel respectively. 
This can be explained by an arrangement of the H atoms in the tetrahedral sites surrounding 
the Sc atoms, thus shielding the Fel sublattice from the Sc sublattice, as already suggested 
by Smit et al.^, meaning a reduction in terms of contact between Sc and Fel. Regarding 
the already mentioned role of the Sc-Fel interaction (see section IV. A. 3) in the onset of 
the magnetic moment on scandium, this reduction is concomitant with a more pronounced 
magnetic moment for Fel and a less large one for Sc -in agreement with the computations- 
. The magnetic moment values given in table [III for Fel and Fe2 are 2.486 and 1.967 eV 
respectively. This inversion for the order of magnitudes, was already observed and explained 
for the expanded hydrogen free ScFe2 (see section IV.B.l). 
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B. Hyperfine field 

The Fermi contact term of the hyperfine field calculated values for the dihydride are such 
as: if|,°*f'(Fel) = -203 kGauss and Hp}?^{Fe2) = -232 kGauss. The i7|?(5^ computed 
values are such as —272 and —214 kGauss for Fel and Fe2 respectively. Experimental 
findings in RefJ of —239 and —300 kGauss correspond to Heff. These values were given 
without being assigned to the two crystallographic sites (6h and 2a) for iron. Since the 
Fel sublattice is shielded by hydrogen from the Sc sublattice, the contact between Fel 
atoms and the neighboring atomic species is reduced, consequently it can be assumed that 
H'^^'^ contribution is the one measured by the experiment for Fel. It follows that the 
largest experimental value of -300 kGauss for -f^e// can be attributed to Fel for which 
the calculations have given a H'^^^ equal to —272 kGauss. On the contrary, Fe2 atoms 
interactions with neighboring species are more important, hence the measured hyperfine 
field corresponds to the sum of both core and valence contributions of Hpc- Then the 
experimental value of —239 kGauss is assumed for Fe2, in agreement with the computed 
value of —232 kGauss for H^p^\Fe2). On the other hand, the orders of magnitudes for Hpc 
(see table [IT]) show that the highest value of Hpc corresponds to the larger nipe magnitude. 
Considering the core part for Fel on the one hand and the total contribution for Hpc of Fe2 
on the other hand, we can propose that this tendency is preserved upon hydriding. 

C. Projected density of states PDOS 

A visual inspection of the spin projected PDOS of ScFe2H2 represented in Fig.[6]shows two 
small energetic intervals ranging from —8 to —5.5 eV and from —11 to —10 eV respectively. 
The similar skylines between the different atomic species featuring in these two energy 
regions describe the hybridization of the metallic species with hydrogen. The sharper and 
narrower nature of these PDOS peaks, compared to those of ScFe2 [Fig. Hl^a)] and the 
expanded ScFe2 [Fig. 111(b)], point to a larger localization of the states in the dihydride 
system. This is due to hydrogen intake and not to volume expansion. 
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FIG. 6: Site and spin projected PDOS of ScFe2H2. For the sake of clear presentation, hydrogen 
contribution was multiphed by 10. 

D. Chemical bonding 

To find an explanation for the magnetic behavior, chemical bonding plots with ECOV 
criterion (shown in fig. [7]) for hydrogen interactions with the metal species from which the 
lattice is built are analyzed. In as far as Hpc originates mainly from the spin polarization 
of valence states via interactions with the magnetic neighbors, the inversion for the order of 
magnitudes mentioned earlier for Hpc is expained by the fact that both hydrogen insertion 
and shielding (see section V.A) reduce Fel interaction with neighbors. Fel atoms become 
further isolated upon hydriding, and tend towards magnetic behaviors of the same kind as 
in Fe metal. In this so-called weak ferromagnetism the majority spin subband as well as the 
minority spin subband are partially depleted. ^^ This is plotted in Fig. [7| where the Fel-H 
interaction is anti-bonding through all the VB. This anti-bonding behavior confirms this 
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FIG. 7: Chemical bonding: non-magnetic ECOV within ScFe2H2 for the atomic pair interactions 
Fel-H, Fe2-H and Sc-H. 

reduction in contact between Fel and the other atomic species. 

VI. ScFe2 AND ANISOTROPY CHANGES 

For the sake of addressing anisotropy effects, an additional expanded hydrogen free ScFe2 
model was calculated with the experimental value of the c/a ratio of the alloy system.^ The 
computed magnitudes and signs for the magnetic moments as well as for Hpc are reported 
in table [TTl An analysis of these results demonstrates that the change of the c/a ratio does 
not affect the general trends of the magnetic behavior. Moreover, the already observed 
inversion (upon volume expansion) for the order of magnitudes for both magnetic moments 
and Hpc is conserved. It can be then concluded that the anisotropy changes are of negligible 
influence on the magnetic behavior of ScFe2H2, i.e., with respect to both volume expansion 
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and hydrogen insertion effects. 

VII. CONCLUSIONS 

In this work local spin density functional (LSDF) investigations of the hydrogen insertion 
effects on the magnetism and bonding within the C14 ScFe2 laves phase have been under- 
taken. In order to address these features, we performed ab initio all electrons computations 
of the electronic band structure and of the bonding properties for ScFe2 and its dihydride 
ScFe2H2 as well as for the expanded hydrogen free ScFe2 system at the dihydride volume. 
Contrary to former studies^'^'^I which described the magnetic behavior of ScFe2 by means of 
a rigid-band shift, our results point to a "covalent magnetism" -like behavior. The analysis of 
the electronic structures and of the chemical bonding properties using the covalent bond en- 
ergy (ECOV) criterion leads to suggest that the volume expansion increases the magnitudes 
for both magnetic moments and Hpc- As for the chemical effect of hydriding, it further 
enhances the magnitude of the total magnetic moment, while the magnetic moments of the 
different atomic species show different behaviors. Furthermore, Hpc values exhibit orders of 
magnitudes which are proportional to those of the magnetic moments of Fel and Fe2 within 
ScFe2, i.e., both m(Fe2) and Hp^\Fe2) are larger than m(Fel) and Hp^\Fel). Volume 
expansion introduces an inversion for this order of magnitudes but respects the proportion- 
ality. However, upon hydriding, one should consider the core and valence contributions to 
Hpc- As a matter of fact, inserted H atoms shield the Fel sublattice, thus reducing contact 
with other sublattices. This leads to experimentally measured values of the hyperfine field 
for Fel concomitant with the calculated H'^^'^ which is proportional to the magnetic moment 
of the onsite rf-electrons. On the other hand, in as far as the hydriding process introduces 
an anisotropic expansion to the ScFe2 unit cell, further computations were performed for an 
expanded hydrogen free ScFe2 model at the experimental c/a ratio of ScFe2 alloy system. 
The analysis of the results of these calculations helped showing that the anisotropy changes 
have no significant effects on the general magnetic behavior. 
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